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In this talk Q] we report work on the matching of the next-to-leading logarithmic ap- 
proximation (NLLA) onto the fixed next-to-next-to-leading order (NNLO) calculations 
for event shape distributions in electron-positron annihilation. Furthermore we present 
preliminary results on the determination of the strong coupling constant obtained using 
NLLA+NNLO predictions and ALEPH data. 

1 Introduction 

The reaction of e + e _ annihilation into three jets has played historically a very prominent role 
for phenomenology. It permitted for example the discovery of the gluon and the measurement 
of its properties and allows also a precise determination of the strong coupling constant a s , 
since the deviation from two-jet configurations is proportional to it. Not only jet rates, but 
also the shape of the single events can be studied in a systematic fashion. The so-called 
event shape observables became very popular mainly because they are well suited both for 
experimental measurement and for theoretical description since many of them are infrared 
and collinear safe. The main idea behind event shape variables is to parameterize the energy- 
momentum flow of an event, such that one can smoothly describe its shape passing from 
pencil-like two-jet configurations, which are a limiting case in event shapes, up to multijet 
final states. At LEP a set of six different event shape observables were measured in great 
detail: thrust T (which is substituted here by r = 1 — T), heavy jet mass p, wide and 
total jet broadening B\y and Bt, C-parameter and two-to-three-jet transition parameter in 
the Durham algorithm y3. The definitions of these variables, which we denote collectively 
as y in the following, are summarized in [3J. The two-jet limit of each variable is y — > 0. 
Until very recently the theoretical state-of-the-art description of event shape distributions 
was based on the matching of the NLLA [3] onto the NLO [H IS] calculation. Using these 
predictions the largest contribution to the error in the determination of the strong coupling 
constant came from theoretical scale uncertainties. Recently the NNLO corrections became 
available. Using this new results we computed the matching of the resummed NLLA onto 
the fixed order NNLO. 

2 Fixed order and resummed calculations 

At NNLO the integrated fixed order differential cross section 

Chad Jo «X 

is given by 

R (y, Q, M) = 1 + «s (/") A (y) + a 2 s (/x) B (y, Xp) + a 3 s (/x) C (y, x^) , 
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where a s = a s /(27r) and x M = /VQ- 

Table [T] shows the relevant contributions for the 



computation of the three coefficient functions A., 
B and C. The careful subtraction of real and vir- 
tual divergences is done using the antenna for- 
malism and implemented in a numerical integra- 
tion program. Recently an inconsistency in the 
treatment of large-angle soft radiation was dis- 
covered [SJ. This was corrected (erratum to [5]) 
and it results in numerically minor changes to 
the NNLO coefficients in the kinematical region 
of phenomenological interest here. The correc- 
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tions turn out to be significant only in the deep Table 1: Contributions order by order, 
two-jet region, e.g. (1 — T) < 0.05 (figure [TJ. 

Approaching the two-jet region the infrared logarithms in the coefficient functions become 




10 2 10 10 2 10 



1-T 1-T 



Figure 1: Comparison between old and corrected distributions for r. In the fixed order 
distribution (left) a small difference is visible in the far infrared region, in the matched 
distribution (right) the curves are equal since the resummation becomes dominant in the 
infrared region. 

large spoiling the convergence of the series expansion. The main contribution in this case 
comes from the highest powers of the logarithms which have to be resummed to all orders. 
For suitable observables resummation leads to exponentiation. At NLLA the resummed 
expression is given by 

R(y,Q,(J>) = (l + Cia.) e (Lai( as L)+ a2 ( as L)) ^ 

where the function g\ (a s L) contains all leading-logarithms (LL), g 2 (a s L) all next-to- 
leading-logarithms (NLL) and [i = Q is used. Terms beyond NLL have been consistently 
omitted. The resummation functions gi(a s L) and g2{a, s L) can be expanded as power series 
in a s L: 

L 9l {a s L) = G 12 a s L 2 + G 23 a 2 s L 3 + G 34 a 3 s L 4 + ... (LL) , 
g 2 {a s L) = Gnd s L + G 2 2(i 2 s L 2 + G 33 a 3 s L 3 + . . . (NLL). (1) 

Table [2] shows the logarithmic terms present up to the third order in perturbation theory. 
At the fixed order level the LL are terms of the form a"L™ +1 , the NLL those which go like 
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Table 2: Powers of the logarithms present at different orders in perturbation theory. The 
color highlights the different orders in resummation: LL (red) and NLL (blue). The terms in 
green are contained in the LL and NLL contributions and exponentiate trivially with them. 



a"L n , and so on. Notice that this can be read off the expansion (JTJ) of the exponentiated 
resummation functions. 

Closed analytic forms for the functions g\{a s L) and g 2 (a s L) are available for r and p [9], 
Bw and Bt PZH E] , C [T3] and Y 3 [13] , and are collected in the appendix of [H] . Recently 
also g 3 (a s L) and 54 (a s L) were computed for r using effective field theory methods [T5] , 

3 Matching of fixed order and resummed calculations 

To obtain a reliable description of the event shape distributions over a wide range in y, 
it is mandatory to combine fixed order and resummed predictions. The two predictions 
have to be matched in a way that avoids the double counting of terms present in both. 
A number of different matching procedures have been proposed in the literature, see for 
example [2] for a review. In the so-called i?-matching scheme, the two expressions for R (y) 
are matched. We computed the matching in the so-called In J?-matching [3] since in this 
particular scheme, all matching coefficients can be extracted analytically from the resummed 
calculation. The In R- matching at NLO is described in detail in [3]. In the In R- matching 
scheme, the NLLA+NNLO expression is 

ln(R(y,a s )) = L 9l {a s L) + g 2 (a s L) + a s (A{y) - G n L - G 12 L 2 ) + 
+ a 2 s (b (y) - \a 2 (y) - G 22 L 2 - G 2Z U 

C(y)-A(y)B(y) + ±A 3 (y)-G 33 L 3 -G 3i L^ . (2) 

The matching coefficients appearing in this expression can be obtained from |T]) and arc listed 
in |16j . To ensure the vanishing of the matched expression at the kinematical boundary y max 
a further shift of the logarithm is made 0. 

The renormalisation scale dependence of ([2]) is given by making the following replace- 
ments: 

a s -> a s (/i) , 
B{y) - B(y,fi) = 2p liix ll A(y) + B(y) , 

C(y) C(y,n) = (2f3olnx^ 2 A(y) + 2lTix^ [2 foB (y) + 2 fa A (y)] + C (y) , 

g 2 (a s L) -> g 2 (a s L, [i 2 ) = g 2 (a s L) + — (a s L) 2 g[ (a s L) Inx^ , 

G 22 G 22 (/*) = G 22 + 2[3 G 12 Inx^ , 

G33 ~* G33 (A*) = G33 + 4/3 G 23 lna;^ . 
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In the above, g[ denotes the derivative of <?i with respect to its argument. The LO coefficient 
A and the LL resummation function gi, as well as the matching coefficients GV;+i remain 
independent on \x. 



4 Discussion of the matched distribution 

For the resulting plots of the matched distributions we refer to [16] . The most striking 
observation is that the difference between NLLA+NNLO and NNLO is largely restricted 
to the two-jet region, while NLLA+NLO and NLO differ in normalisation throughout the 
full kinematical range. This behavior may serve as a first indication for the numerical 
smallness of corrections beyond NNLO in the three-jet region. In the approach to the two- 
jet region, the NLLA+NLO and NLLA+NNLO predictions agree by construction, since 
the matching suppresses any fixed order terms. Although not so visible on these plots, 
the difference between NLLA+NNLO and NLLA+NLO is only moderate in the three-jet 
region. The renormalisation scale uncertainty in the three-jet region is reduced by 20-40% 
between NLLA+NLO and NLLA+NNLO. This effect is due to the smaller renormalization 
scale dependence of the NNLO contributions. It is also important to observe that the scale 
dependence remains the same and is larger in the two-jet region, because the resummed 
calculations at NLLA take into account only the one- loop running of the coupling constant. 
This has important consequences in the determination of a s and we will comment more on 
this in the next section. 

The description of the hadron-level data improves between parton-level NLLA+NLO 
and parton-level NLLA+NNLO, especially in the three-jet region. The behavior in the two- 
jet region is described better by the resummed predictions than by the fixed order NNLO, 
although the agreement is far from perfect. This discrepancy can in part be attributed to 
missing higher order logarithmic corrections and in part to non-perturbative corrections, 
which become large in the approach to the two-jet limit. 

The right plot in figure [T] shows that the inconsistency in the treatment of the large-angle 
soft radiation does not affect the matched prediction since the infrared region is dominated 
by the resummation. 



5 Determination of the strong coupling constant 

After the extraction of a s using only the NNLO distributions and the experimental data of 
ALEPH [14j . a new extraction of a s using the new matched results was performed using 
JADE data [18] . The improvement in the error coming from the inclusion of resummed 
calculation is not as drammatic as passing from NLO to NLLA+NLO calculations. As 
already anticipated, this is due to the fact that the NNLO coefficients compensate the two- 
loop renormalization scale variation, whereas the NLLA part only compensates the one-loop 
variation. A more natural way of matching would be the consider NNLLA and NNLO, but 
the NNLLA function (73 is by now only known for r. A new determination of a s using 
ALEPH data is in progress. The analysis will follow the lines of the previous determination 
using pure NNLO predictions with a few improvements. 
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6 Outlook 



The matching of NLLA and NNLO has improved the theoretical prediction of event shape 
distributions, but further improvement is possible by including the NNLL corrections into 
the calculations. These corrections are known only for r, where higher order logarithmic 
corrections have been computed [T5] using soft-collinear effective theory (SCET). From these 
calculations one can extract the functions 173 (a s L) and 54 (a s L). The next step towards the 
further improvement in the extraction of a s from event shape distributions could be to 
compute them for all six observables mentioned here. As shown in |15| the subleading 
logarithmic corrections can also account for about half of the discrepancy between parton- 
level theoretical predictions and hadron-level experimental data. 

Improvements can also come from non-perturbative corrections. A very recent non- 
perturbative study for r using a low-scale effective coupling [TT] shows that non-perturbative 
l/Q power corrections cause a shift in the distributions, which can account for an important 
part of the difference between parton-lcvcl distributions and hadron-level experimental data 
discussed in the previous section. 
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